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ABSTRACT 


Meroe mmouemroreaerermining flight director laws for the 
longitudinal control of a V/STOL aircraft in landing 
approach is evaluated. The method is based on the application 
of an optimal control model for the human pilot. The vehicle 
studied was the UH-1H helicopter at three approach ground- 
speeds: 60 knots, 40 knots, and 20 knots. The two pilot 
Supeues were Jonesitudinal cyclic and collective. In the 
analysis, ten pilot "transfer functions" which relate the 
two control variables to the five displayed and perceived 
quantities were obtained, These transfer functions were 


PaeoeUtLiiazed to obtain the respective flight director laws. 
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ieee NT ROD eTION 


VlOleaAtrCreamemiave recelved considerable attention for 
paramilitary ana Clvyilian use, due to their unique ability 
fmemoperavce In Areas not accessible to conventional aircraft. 
The rotary-wing aircraft, in particular, has shown such 
increased diversity that display systems which will allow 
VTOL operations in zero/zero conditions are being actively 
sought. Well-designed flight directors are the heart of 
such systems, 

A flight director is a display system which provides 
control commands to the pilot to enable him to complete 
demanding flight tasks with relative ease and precision. 
Recently, Levison /Ref. 1/ has proposed an sane aa 
wechnique for flight director design based upon the optimal 
pilot model pioneered by Kleinman, Baron, and Levison /Ref. ae 
Pemogitfied form of this technique is utilized in this thesis 
to obtain the longitudinal flight director laws for the UH-1H 
helicopter in three approach conditions in the presence of 
random vertical and horizontal atmospheric turbulence. Five 
mentele motion quantities, which are normally directly 
Slepravyed or perceived by the pilot, are blended to drive two 
ioolaye symbols, the cyclic director and the collective 
director. The design was undertaken at three different 
groundspeeds: 60 knots, 40 knots, and 20 knots. The pilote- 
vehicle system for the optimal modeling procedure is shown 


ereapumeit ly an Figure 1, 
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II. METHOD OF ANALYSIS 


A. THE MODELING HYPOTHESIS 


Subject to his inherent limitations, the well-trained, 


well motivated pilot behaves in an optimal manner. The 


pilot's control characteristics can be modeled by the solution 


of an optimal linear control problem and an optimal estimation 


problem with certain "modifications." 


B. MODIFICATIONS FOR PILOT MODELING 


pMeoube tinesdelay 15 included in each of the pilot's 
COnproet OULDULS . 


Each output neuromuscular system is modeled as a 
fiest order lag. 


Each observed variable is assumed to contain pilot 
induced additive white noise which scales with the 
variance of the observed variable. Also, each 
control output is assumed to contain pilot induced 
additive white noise which scales with the variance 
of the control motion. 


If a variable is displayed explicitly, the pilot 
also perceives the first derivative of the variable 
but no higher derivatives. The first derivative is 
also noise contaminated. 


Micmmaeeoxn Of Performance for the optimization 
procedure is chosen subjectively to mirror what the 
display system designer believes to be the task and 
control objectives perceived by the pilot. 


C. LONGITUDINAL HELICOPTER EQUATIONS OF MOTION 


The longitudinal helicopter equations of motion are 


Shown in state variable form on the next page. The 


assumptions used in the derivation of these equations follow 


11. 
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Meerneavenmete as idealized aS a rigid airframe to which 
is avtached a rotor. 


2. The rotor is described by the tip path plane whose 
orientation determined propulsive and aerodynamic 
forces and moments. 

3. No rotor degrees of freedom are considered other than 
control inputs which serve to describe instantaneous 
tip path plane orientation. 


4. All coupling between longitudinal and lateral 
motion is ignored. 


5. There is linearized small perturbation motion about 
the horizontal reference flight path. 
D, ATMOSPHERIC TURBULENCE 
The power spectral densities for the atmospheric turbulence 


are those suggested by Hart, Adkins, and Lacau [{Ref. Cre 


20 wlny 1 
Pos W 2) Ng U Inyo , 2 oD 
gs O i (>) 
O 
2 
(w) = Ae Nes aL 
e"g U5 (nee ce 
U 
O 
with 
2U. Un 
a 1 oe B, a fone (8) 
W W 
2U . Un 
NT oe BZ, °= a (9) 
u U 


Using the concept of a white noise excited shaping filter, 


the following state equations can be developed. 


Wo “BW, + a.V (10) 


ea | eae 
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i, = -B,u, + o,v, | (11) 
Equations (10) and (11) now augment equations (1) - (5), 


E, NEUROMUSCULAR AND TIME DELAY EQUATIONS 

As noted in section II.B., a pure time delay and first 
ordger lag are included in the pilot's output. The delay is 
approximated by a second-order Pade' function. The quality 
of first and second order Pade' approximations are indicated 
in Figure 2 for the frequency interval of interest in pilot 
modeling, wt < 5 RAD. For ease of development, a portion of 


Figure 1 is shown here. 
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COMPARISONS 














‘ See 16 ‘ 8 16 
Ud t+ 7 Uqt Te Yoda = Us a a (14) 
Now let 
a> Ug aan Gils) 
8 16 
ay = ti ar u, (16) 
° 16 
we —— v, (17) 
Je 16 8 128 
d, = d, a d, ate a , (18) 
il 1 i) Oa 
6, =-a- d, - = 6, + wm Uu, + (19) 
B Ty 1 Ty B Ty 1 We 
Using the same technique to develop ae cine, and on as was 
used to develop a dg and ae one obtains: 
a 16 
d, = d, = T us (20) 
e —=-16 8 128 
il f 1 oe 
6.=—d,-=—6,+ =u, + (22) 
6 =i 3 ies C sae 2 ie 
Pemevions (1) - (5), (10) - (11), and (15) - (22) now 


constitiute the state description of the pilot-vehicle 


system, 


F, DISPLAYED VARIABLES 
The displayed variables listed below are those quantities 


which are assumed to be displayed or perceived by the pilot, 
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Zi = Displayed Groundspeed Deviation = K.u 


Dp, : 
Z = Perceived Time Rate of Change of Glideslope 
P» Deviation = ew + Un OQ ) 
2 = Perceived Pitch Rate = K q 
3 3 
Zi = Displayed Pitch Angle Deviation = K,0 
4 
2. = Displayed Glideslope Deviation = Koh 
5 


The Ky are display gains. For example: 


| radians subtended at the pilot's eye by display 
3 rad/sec pitch 


element motion 
rate 


Table 1 lists the display gains for this study. 

meccording to modification (4) in section I1,B., u 
Should also be a perceived variable. However, since all 
perceived variables will be used in generating the director 
laws to be discussed, they must be measurable. If u repre- 
sents airspeed, u will be difficult to sense. For this 
reason wu was not considered a perceived variable (despite 
the fact the u represented groundspeed in this particular 
analysis), 

It should be emphasized that neglecting u imposes little 
constraint on the model's validity since this variable 
1s associated with the phugoid and will be quite small 


Eoroughent the approach. 
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TABLE I. Display Gains 


G. OBSERVATION AND MOTOR NOISE 

At present it is not possible to Separate experimentally 
the various sources of pilot "remnant." In this paper, as 
is done in reference 2, observation noise and motor noise 
were taken to be representations of remnant. Observation 
noise was taken to be pilot injected noise in observing the 
displayed quantities, and motor noise was considered to be 
random error in control movement execution. Both were 
considered independent, Gaussian, white noise processes, 

aie pfrOL related noise Levels were set at values 
considerably larger than those found in documented laboratory 
experiments (e.g. Ref. 2) so that the design would be less 
sensitive to pilot noise and more forgiving of actual non- 


optimal pilot behavior. Consequently, numerical values for 
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the observation noise and motor noise were chosen as p = 0,1, 


and o' = ,Ol1. 


H,. SYSTEM EQUATIONS AND MATRIX NOTATIONS 
The system equations listed below define the optimal state- 
feedback controller and estimator problem. A thorough develop- 


ment of these equations can be found in reference 4. 


—(t mew oe(t) + Bu (t) + y we(t) (285 
y(t) = C x (t) | (24) 
aye= Wx (t) + v (t) = 2. (t) + v (t) C25) 
A: 
E | wot dw (T + t)]= ac (26) 
E [vctovs (a + «| = G O(T) C2) 
di 
pe nea f [xe ay ce eure) Ru ce] at (28) 
O 
where 


ees an nx nm plant matrix 

Pee) cea x 1 state vector 

Peels alex p Control matrix 

tee) 1S 2 px 1 control vector 

ilocos mn xX t disturbance matrix 

w(t) is a t x 1 disturbance vector 

Vieeweicea q x 1 output vector 

Ceicmand < 1 OUTpUt Matrix 

Hee Cty tS a vector of linearly uncorrelated, zero 


mean white noise signals with Gaussian amplitude probability 
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fastribution functions, The elements of w(t) are assumed to 
be sample functions from n random processes which are each 
ergodic and are jointly ergodic. The covariance matrix for 


w(t) is 
i | w(t) wi(t + t)] = Ewoqr) (29) 


where 6(tT) is the unit impulse function. 


The measured quantities on sensor Signals are 


z(t) = H w(t) + v (t) | (30) 


where 
Aga 1oea vex I measurement vector 
H is a u X n measurement matrix 
G2 sts tee Measurement noise vector 
The elements of v(t) are assumed to be sample functions from 


p random processes each of which are ergodic and jointly 


eugodr1c, The covariance matrix for v(t) is 


E | v(t) vi(t + =) | = 6 él) (31) 


The system is assumed to be completely controllable and 
completely observable. It is desired to find the control 


mietion u(t) which minimizes the quadratic scalar index of 


performance 
i 
J = ae = f | v7Ct) Ory CG) Fr u'(t) R u(t) | at (32) 
O 
where 
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Q is a q X q Symmetric output cost weighting matrix and 
at least positive semidefinite 
Ris ap xX p Symmetric control cost weighting matrix and 


positive definite 
The solution to the linear quadratic Gaussian control 
problem can be outlined as follows: 


a.) The optimization problem can, by the called Separation 
Theorem, be broken up into two separate problems, an 
optimal control problem and an optimal estimation 
Or fiiberins problem. 


b.) The optimal estimation or filtering problem generates 
an optimal estimate, X(t) of the state x(t). This 
estimate is optimal in the sense that 


To 
lim 1 X(t) K(t) at 
7S 
O 


'T'->00 


is minimized, where x(t) is the estimation error defined as 
AG) Secs) = Ge) | C3) 


The optimal estimator (or Kalman filter) has the form 


| >> 


A(t) = A R(t) + Bu(t) + K[2(t) - H RCL) (34) 


The estimator gains are given by 
K= PH G™ (35) 


where iP TSypme error covariance matrix 
a Pb x 
E [ R(t )s Cte =) | = P 6(t) | (36) 


Peewee post live detinjte Solution fo the steady — stace 
Pelver mem Nn TCeCall Cqulabion 
AP+PA'+yFy -PH GHP =0 (37) 
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ec.) The optimal control problem generates an optimal 
Cito wlave wt jewiiehn 1s a linear function of the 
estimated state 


u(t) = - L &(t) (38) 


maere L is a p X n optimal controller gain matrix. The gain 
ere x i 1S mdentecal to tne one obtained by solving the 
optimal control problem with no system disturbance, exact 
state information, and the index of performance given by 
s= f [sty @ y(t) + ute) R uty] at (39) 
O 


mooeconuroller@eain matrix L is given by 


L=R" BS (40) 


wnere S 1s the positive definite solution to the steady-state 
control matrix Riccati equation 


ewes Oe s BR BS = 0 (41) 


pecan be Snown that the state covariance matrix 


B [x(t) xi(e + t)] = (B+ MD 6(t) (42) 


Gee to umers@llueion ot the filer matrix Rieccati equation 


and M is the positive definite solution to 


(A- BL) M+M(A-BL)' +KGK’" = 0 (43) 


In addition to the solutions outlined above, it can be 


Shown that the transfer matrix relating the Laplace transform 
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Gf the optimal control law u(t) route wuaplace transtorm of 


vee medsurement vector Z(t) (with v (t) = 0) is given by 

U (S) = -L (SI-A+BL+KH)™ K Z (8) (44) 
where 

v (s) =f[u (| (45) 


IN 


(s) =f z (ty) (46) 


The state variables for this study were chosen in the 


following order 





om, Gidm GS Cle ee Ome ee 


OF WD DnYdD FB wT 


The following matrix tables were deveioped from equations 
ioe mee ovn( Om ( Ll). andes (17) —s(22)),) The matrices are 
tabpeled in accordance with equations (23) = (27), 

As Table VI indicates, all but one of the elements of 


F and G are dependent upon variances of system variables 
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which are not Known axpriori, Estimates of these variances 
must be made, the solution to the optimal estimation and 
eomerol problem obtained, and the resulting variances used 
in a second iteration. This iterative process continues 
until the equations for the 5 and Vs. in Table VI are 
i 

satisfied. As reference 2 points out, 2*m iterations are 
usually involved, where m is the number of displayed and 
perceived quantities. 

ieee ptinemveecvor y utilized in the index of performance 


is given in Table V as 


The elements of the index of performance weighting 
Maerereces Q and Rowere chosen as the reciprocals of the 
"maximum allowable" deviations of the output and controls. 
mais, when an output or control variable attains these 
subjectively chosen magnitudes, it makes a contribution of 
unity to the integrand of the index of performance. The 
maximum allowable deviations were chosen subjectively and 


varied with the flight conditions as Table V indicates. 
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III, PILOT MODELING EXAMPLE 


This section presents a numerical illustration of the 
pilot modeling technique developed in Section H. The 
longitudinal dynamics of a UH-1H helicopter at an approach 
groundspeed of 60 Knots are utilized. The normalized UH-1H 
longitudinal derivatives, ina Stability axis system, are 
Shown in Table VII, 

A modified form of the Variable Automatic Synthesis 
Program (VASP) /Ref. 5/ was utilized to solve the optimal 
estimation and control problem. The solution tc this 
problem constitutes the pilot-vehicle model. After roughly 
10 iterations the solution for the 60-knot case converged. 
Table VIII shows the root-mean-Square (RMS) performance 
figures. The ten pilot transfer functions relating the cyclic 
and collective control variables to the five displayed and 
perceived quantities were obtained. Figures 3 and 4 are 


Bode plot representations for the (s) and a (s) pilot 
h 


transfer functions respectively. 

Pure gain approximations were then made to each of the 
transfer functions in the frequency range of interest for 
modeling purposes: .1< w < rad/sec. The gains were then 
normalized by dividing by the magnitude of the largest gain 
hor Cach control, The values in Table IX Poniiteds Eor 
example uSing the gains of Table IX the cyclic and collective 


flight director laws become 
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é,(t) = - 8.42 « Hen Gee eee 107 h(t) = ace) 


Peg45 ©(t) 2 4002 « 10° "h(t) (47) 
. 7 a =Sic 
da(t) — cenieemlOor wor) = S.02 * LO’ h(t) « q(t) 

~- ,1358 O0(t) + 1.96 «+ omen) (48) 


where Op and 56 Le weneneyvelice anamcollective director 


Signals respectively, 
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TABLE VII. 


2eS-00S-35/T 00 HEO8TL'O = on 
Ho 
99S-99S-13]/T Onde scece G= = i 
99 
.9S-99S/T €0O ApPP6tl’°O- = 7 
Ho 
oes-09S/T COmayeece 0 = 7; 
99 
oes-00s/T 00 HZE8SZE°O- = Xx 
Ho 
098S-090S/T TO HI6SO8 ‘0 = X 
M 
aT 0'O = “HN 
hs 
oes/T TO wSsoSs8sz'O- = W 
M 
ei =09s/ 1 MO BUe ge = = N 
nN 
4j-00S/T ZGO-aALZCSPL O = A 
b 
oes/4J ZO AOVS9E'O = 7, 
M 
0es/T 00 HOLP68'O0- = Z, 
nN 
oes/T TI0-S88e6r°0- = Z 
b 
oes/4F COmmetLoOoST Oo = Xx 
+ OM 
o0s/T  TO-H6686T°O = X 
n 
Oo5/—  eO-a7CIrs 0- = x 
O 
Sas/ at €O HPETOT‘'O = n 


WHLSAS SIXV ALITIGVLS NI SHYAILVAITYAd TVNIGOLIONOT HI-HN GAZITVAYON UL 
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Performance 
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SENSORY LOW-FREQUENCY 
DIRECTOR VARIABLE GAIN 


Creche 
COLLECTIVE 


TABLE IX, Gains for UH-1H Helicopter Director Laws - 









Airspeed -.000842 sec 








Sink Rate -,.000245 sec 


ft-~sec 



















Pitch Rate 1 
rad 
paeen 1945 Lt 
= : rad 
Height ~ 000402 ag 







Airspeed ~000247 sec 


Sink hate ,00552 sec 


ft-sec 
rad 














Pitch Rate -1 










uae 
rad 
tet 


/OO196 F— 





|231 16} 9 -,1358 





Height 


Velocity 60 Knots 
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IV, RESULTS AND CONCLUSIONS 


The normalized director gains for the 40 and 20-knot 
approach speeds are shown in the following tables. It should 
be noted that all of the Bode magnitude plots for the thirty 
Dr1lot transfer functions (ten each for each of the three 
approach speeds) could be approximated by a fifth-order low 
pass filter. The break frequency varied somewhat between 
O74 and 0.5 rad/sec. Only the director control gain changed 
meom Lunction to fumetions 

in order to) ascertain, the reason for the rather dramatic 
Sign reversals which occurred in some of the larger gains as 
the flight condition changed (e.g. in the collective-to-pitch 
rate and collective-to-pitch gains of Table IX and Table X), 
the author re-ran the 40-Knot case with identical pilot 
Po eemneters  spectiically @ and R matrices, as in the 60-knot 
case. The resultS were Similar to those of Table X. This 
manmea tea thatetne flight condition dictated the gain sign 
ilar lOn saimdmmot the Subjecvive Selection of the Q and R 
matrices in the pilot model. 

The director laws implicit in Tables IX ~- XI must be 
evaluated in piloted simulation before the efficiency of the 


design method outlined in this thesis can be determined, 








SENSORY LOW-FREQUENCY 
DIRECTOR VARIABLE GAIN 


Airspeed ,VO0D3G.SeC 
Sink Rate ~,Q005287 sec 
CYCLIC S1t@n Vane 1,9 Lt-sec 
rad 
Pitch 374 2h 
rad 
ft 
Height -,00175 Ft 
Airspeed .000426 sec 
Sink Rate -~,0044 sec 
COLLECTIVE Pitch Rate 10 Ease’ 
upeln 41962 12_ 
rad 
ft 
Height -,00188 == 


TABLE X. Gains for UH-1H Helicopter Director Laws - 





Velocity 40 Knots. 
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SENSORY LOW-FREQUENCY 
DIRECTOR VARIABLE GAIN 


Airspeed -.0126 sec 
Sink Rate 004448 sec 
CYCLIC Pitch Rate i) SSeS 
rad 
Pitch oa SE 
rad 
, tt 
Height "OOZ99 Tt 
Airspeed ~,01257 sec 
Sink Rate .00435 sec 
COLLECTIVE Essel. ever yj 2tcsee 
rad 
Pitch co 
rad 
Height 002535 rf 


TABLE XI, Gains for UH-1H Helicopter Director Laws « 








Velocity 20 Knots. 
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